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Quark Confinement from Color Confinement 
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We relate quark confinement, as measured by the Polyakov-loop order parameter, to color confine- 
ment, as described by the Kugo-Ojima/Gribov-Zwanziger scenario. We identify a simple criterion 
for quark confinement based on the IR behaviour of ghost and gluon propagators, and compute the 
order-parameter potential from the knowledge of Landau-gauge correlation functions with the aid 
of the functional RG. Our approach predicts the deconfinement transition in quenched QCD to be 
of first order for SU(3) and second order for SU(2) - in agreement with general expectations. As an 
estimate for the critical temperature, we obtain T c ~ 284MeV for SU(3). 
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INTRODUCTION 



Aside of the confinement of quarks, the confinement 
of gluons is a challenging and unsolved problem. Vari- 
ous scenarios predict the confinement mechanism to be 
manifest in the infrared domain of gauge-dependent cor- 
relation functions. In the Kugo-Ojima [l[ and Gribov- 
Zwanzigcr scenarios [2[ (KOGZ) an infrared enhancement 
of the ghost and an infrared suppression of the gluon sig- 
nal confinement. These scenarios have been investigated 
by a variety of non-perturbative field theoretical tools 
such as functional methods [H, 0, Q and lattice gauge 
theory The results provide strong support for these 
scenarios even though the infrared enhancement of the 
ghost is a subject of ongoing debate, for a summary see 
e.g. This paves the way for a comprehensive under- 

standing of the non-perturbative mechanisms of strongly- 
coupled gauge systems. 

A pressing open question is the relation of color 
confinement to quark confinement. Typical quark- 
confinement criteria based on the Wilson-loop or 
Polyakov-loop expectation value [n| in quenched QCD 
have so far remained inaccessible from the pure knowl- 
edge of low-order correlation functions of the gauge sec- 
tor, although evidence for a linearly rising potential be- 
tween static quarks has been collected within certain ap- 
proximation schemes, e.g. fill. [Til. . 

In this letter, we propose a method for computing 
the full Polyakov-loop potential from background-field- 
dependent Green functions. Our approach relates the 
order parameter of quark confinement, the expectation 
value of the Polyakov loop, to the momentum depen- 
dence of gauge-dependent Green functions. This leads 
to a simple confinement criterion in any gauge. The 
method is explicitly applied in the Landau gauge, where 
it relates the KOGZ scenario of gluon confinement to 
quark confinement. We evaluate the effective potential 
of a purely temporal background field configuration Aq, 



being directly related to the Polyakov loop variable, 



L(x) = — trP cxp 



iff 



dx A (x ,x) 



(1) 



where V denotes time ordering, and the group trace is 
taken in the fundamental representation. The negative 
logarithm of the Polyakov loop expectation value relates 
to the free energy of a static fundamental color source. 
Moreover, (L) measures whether center symmetry is re- 



alised by the ensemble under consideration, see e.g. [14 1. 
A center-symmetric confining (disordered) ground state 
ensures (L) = 0, whereas deconfinement (L) ^ signals 
the ordered phase and center-symmetry breaking. 



The order parameter (L[^4o]) is conveniently 
parametrised in the Polyakov gauge: doA = 
with Aq in the Cartan subalgebra. Then, (Aq) is 
sensitive to topological defects related to confinement 
[H|], and also serves as a deconfinement order pa- 
rameter. More specifically, (L[Ao]) is bounded from 
above by L[(^4o)] owing to the Jensen inequality 
L[(A )} = trexp(ig /3(A ))/N C > (L), such that L[(A )} 
is nonzero in the center-broken phase. In the center- 
symmetric phase where the order parameter (LLAo]) 
vanishes, also the observable L[(^4o)] can be shown to be 
strictly zero [l6j]. This establishes both (Aq) as well as 
L[(Ao)] as a deconfinement order parameter. 

In the present work, we compute the effective poten- 
tial for (Ao)Jrom Green functions in the background-field 
formalism [l7| in the Landau-DeWitt gauge by means of 
the functional RG. These Green functions can be deduced 
from that in the Landau gauge, that is at vanishing back- 
ground field. Our construction relates gluon confinement 
encoded in the IR behaviour of Green functions to the 
potential of the order parameter for quark confinement, 
and provides a simple confinement criterion. 



2 



BACKGROUND-FIELD FLOWS 

The effective potential is given by F(L[Ao]) = T/fi, 
where T is the effective action taken at the mean field Ao, 
and £1 is the space-time volume. We evaluate the effective 
action T in the background field approach, where T on 
the one hand depends on the field variable A, being the 
expectation value of the fluctuating quantum field. On 
the other hand, a dependence on an auxiliary background 
field A is introduced by gauge-fixing the fluctuating field 
with respect to the background, 



D^(A)(A-A) fl = 0. 



(2) 



Implementing this gauge condition at vanishing gauge 
parameter constitutes the Landau-DeWitt gauge. With 
the gauge fixing ([2|) , the field dependence of the effective 
action can be summarized as, T = r[$, A] with fluctua- 
tion fields $ = (A — A, C, C) relative to the background. 
The important connection to the standard effective ac- 
tion depending only on A is established through the iden- 
tity T[A] = T[0, A = A], 0. 

In the present study, we identify the background field 
with the Polyakov loop field, A = Ao. For evaluating the 
effective potential V(L[Aq\), it suffices to consider Aq as 
constant, yielding 



V k (L[A }) 



T k [0,A 



(3) 



Wc compute the effective potential non-pcrturbatively by 
means of the functional RG (FRG) for the effective action 



regulator (|5|) the flow of the standard effective action 
r^[yl] = T k [Q, A = A] is also gauge invariant. 

The flow of r^A] can, in principle, be obtained from 
Eq. (g]) by setting $ = and A = A. But, this flow is 
not closed [5^, (22|: the right-hand side of (g]) depends 
on T^' 0) [Q,A] ^ r^ 2) [A], the flow of which cannot be 
extracted from t^r^lyl]. This has been neglected in pre- 
vious non-perturbative applications [25| but turns out to 
be crucial for confinement. Hence the key input, the two- 
point function r£ 2 ' ' ) [0, A] in the background field, has to 
be computed separately. 



EFFECTIVE ACTION FROM LANDAU-GAUGE 
PROPAGATORS 

First, we observe that in the Landau-DeWitt gauge 
the longitudinal components of Green functions decou- 
ple from the transversal dynamics, which further reduces 
the truncation error, for a detailed discussion see 0. 
Moreover, rjj 2,0 ^ [0, 0] (p 2 ) corresponds to the propagator 
in the Landau gauge, since the background field gauge 
with gauge condition ^ reduces to the Landau gauge 
for vanishing background field. The Landau-gauge prop- 
agator has been computed within functional methods, 
[3, H 0|, as well as within lattice gauge theor y B| ; for 
reviews and further literature, see [a. 1 14, Il9l. l2dL l27l|. 

Recalling the results for Landau-gauge propagators, 
the gluon propagator can be displayed as 



18[, for reviews see g^l^. The flow equation for A] r^ 2 ' 0) [0, 0}(p 2 ) = p 2 Z A (p 2 )P T (p)H + p 



in the background-field approach reads 



, Zl(p 2 ) 



M (6) 



d k T k [$,A\ 



iTr- 



1 



where 1 T 



(n,m) 



r ( *' 0) {<f>,A}+R k 



;T k |2fJ, |21|, |2S 



-d k R k 



(4) 



where n L ^„(p) = p^/p 2 , P T = 1 - P L , H afc = S a b, and 
£ denotes the gauge parameter. For the ghost, we have 



k ~ sWJa^ 1 fc ^ ^ ^J- The regulator 
function R k implements an IR regularisation at p 2 ~ k 2 , 
and the trace Tr sums over momenta, internal indices 
and species of fields. The flow ((4]) interpolates between 
the classical action in the UV and the quantum effective 
action T = T k=Q in the IR. For $ = 0, Eq. (j4|) entails 
the flow of Tfefyl] = r/jO, A = A], and as a specifically 
interesting case, that of T4(L[ylo]) = rfc[Ao]/£l. 

Background-field flows have been applied successfully 
to non-perturbative analyses of chiral properties in full 
QCD [23} , including quantitative estimates of the critical 
temperature of the chiral transition from first principles. 

The flow (|1]) is solved utilising optimisation ideas pol . 
24 | that minimise the truncation error. Here, we use a 
specific optimised regulator [2C| , 



rp°\o,o)(p 2 ) = p 2 Zc(p 2 )v. 



(7) 



The longitudinal dressing function obeys = 1 + (D(£) 
and hence drops out of all diagrams beyond one loop in 
the Landau gauge £ = 0. The dressing functions Za,c 
encode the nontrivial behavior of the full propagators. 
In the deep infrared, they exhibit the leading momentum 
behaviour 



Z A (p 2 - 0) ~ (p 2 ) 



2\K A 



Z c {p 2 -+ 0) ~ (p 2 ) 



2\K C 



(8) 



i? op t = (fc 2 - if } [0,A])9(k 2 



ifioMD, 



(5) 



supplemented by fc-dependent fields $ such that 
r^ 2) [0,i] = r^ 2 2 [0,A] for r^ 2 2 [o,A] > k 2 . With the 



In the last years it has become clear that Landau gauge 
Yang-Mills admits a one-parameter family of infrared so- 
lutions consistent with renormalisation group invariancc 
Despite some formal progress the full understand- 
ing of the underlying structure is a subject of current 
research. Technically, the parameter corresponds to an 
infrared boundary condition, the value of Zc{0), and is 
also relates to Za{p 2 — > 0) [Sj. This fact is reflected in 
recent lattice solutions [29j and indications thereof have 
also been seen in the strong coupling limit (30j . For 
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FIG. 1: Momentum dependence of the gluon (left panel) and ghost (right panel) 2-point functions at vanishing temperature. 
We show the FRG results from Ref. [i| (black solid line) and from lattice simulations from Ref. [f| (red points). 



Zc(p 2 — > 0) — ^ it can be shown that there is a unique 
scaling solution, (3ll . l32l |. Then the two exponents are 
related and obey the sum rule 



= k a + 2k c + 



(9) 



in d dimensional spacetime [J, |28|, [3l| . Possible solutions 



are bound to lie in the range kc £ [1/2 , 1], see |28(. For 
the truncation used in most DSE and FRG computation, 
we are led to 



kc = 0.595... and ka = —1.19. 



(10) 



being the value for the optimised regulator [f|. The reg- 
ulator dependence in FRG computations leads to a range 
of kc € [0.539 , 0.595], see [Ej]; for a specific flow, see [33j | . 
These results entail the KOGZ confinement scenario: the 
gluon is infrared screened, whereas the ghost is infrared 
enhanced with Kc > 1/2. 

In turn it can be shown that for non- vanishing Zc(0) 
the gluon propagator tends to a constant in the infrared, 
^ZAf)-* m 2 , for related work see e.g. @, [H, [H, 
371 . l38l . l39j . Note that the gluon propagator then does 
not correspond to the propagator of a massive physical 
particle. Instead, we observe clear indications for posi- 
tivity violation in the numerical solutions for the gluon 
propagator related to gluon confinement, @,41|. Still the 
gluon decouples from the dynamics as docs a massive par- 
ticle, hence the name decoupling solution. The value of 
Zc seems to be bounded by its perturbative value from 
above, and the gluon mass parameter is bounded from be- 
low |8| . The qualitative infrared behaviour is then given 
by the infrared exponents 



K A 



1 



and Kc = • 



(11) 



We emphasise that even though the infrared exponents 
for the scaling solution (flQ|) and the decoupling solution 
(fTTjl arc rather different, the propagators do only differ 
in the deep infrared. It has been suggested in [8j that 



the infrared boundary condition is directly related to the 
global part of the gauge fixing, and hence to different 
resolutions of the Gribov problem. Indeed in [29[ the 
infrared boundary condition has been implemented di- 
rectly as a global completion of the gauge fixing. Note 
also, that for Landau gauge Yang-Mills with standard 
local BRST invariance the requirement of global BRST 
singles out the scaling solution. The existence of such 
a formulation on the lattice has been shown recently in 
j42l ] . In summary the results are affirmative for the above 
interpretation and are supported by results in the strong- 
coupling limit [30l | for different implementations of lattice 
Landau gauge. 

In turn, it has been also shown in a series of works 
that an infrared condition also is present in Landau g au g e 



37, 3 



Yang-Mills with the horizon function, e.g. [35j,[36 
The latter introduces an explicit (or soft) breaking of 
BRST invariance as it restricts the functional integral 
to the first Gribov region. Still this does not fix global 
gauge degrees of freedom as also the first Gribov region 
contains infinite many gauge copies. The possibility of a 
scaling solution in this framework hints at the validity of 
Zwanziger proposal: full BRST invariance is recovered in 
the thermodynamic limit if the path integral is restricted 
to the fundamental modular domain with only one gauge 
copy. 

In summary a consistent picture has emerged with 
nicely relates all current results. The confirmation of 
this picture certainly would provide further insight to 
the confinement mechanism. For the present work, we 
simply note that the scaling solution is singled out by 
global BRST invariance which allows the construction 
of a physical Hilbert space from gauge fixed correlation 
functions. Nonetheless, the whole one-parameter family 
provides consistent gauge-fixed correlation functions of 
Yang-Mills theory and physical observables should be in- 
sensitive to the parameter choice. In the present work, 
we can test this statement. 

We proceed by extending the Landau-gauge propa- 
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gator to that in a given background A. The Landau- 
gauge two-point function T^'°\o,0]{p 2 ) is, apart from 
its Lorentz structure provided by the projection opera- 
tors Pt/l(p)> a function of only the momentum squared 
p 2 , cf. Eq. ©. At vanishing temperature, the back- 
ground field propag ator if' 0) [0,A] can be related to the 
Landau-gauge propagator in a unique fashion owing to 
gauge covariance, 

(Tj?' 0) [0,A])# = (T^°\0M-D 2 ))f v + F^ v f a (D), 

(12) 

with non-singular /(0) in order to ensure the proper limit 
of a vanishing background. The projection operators 
P T / L implicitly contained in r^ 2,0 ' [0, 0] (— D 2 ) generalize 
to projectors on transversal and longitudinal spaces re- 
spectively with respect to the covariant momentum D, 
Pt/l = Pt/l(D). The / terms cannot be obtained from 
the Landau-gauge propagator, but arc related to higher 
Green functions in Landau-DcWitt gauge. However, for- 
tunately they do not play a role for our purpose. 

At finite temperature, the Polyakov loop L is a further 
invariant, and the 00 component of the gluon two-point 
function (|12[) receives further contributions proportional 
to derivatives of L. For constant fields we arrive at 



(Tp 0) [0,Ao]) a ; v = (rriO,0](-Z^))^+L-tcrms, (13) 

as the / term in JT2]) vanishes: F(Aq) = 0. In this letter, 
we take only the explicit T dependence due to Matsubara 
frequencies into account and drop any implicit T depen- 
dence: first, this amounts to dropping the L contribution 
in (1131). This term is related to the second derivative of 



.(2,0) 



the effective potential vjf^ via Nielsen identities 2(J 2 



and can indeed be estimated by V^ 2 \ Its influence on V 
the confinement-dcconfincmcnt phase transition temper- 
ature is paramctricaUy suppressed, and can be neglected 
for a first estimate of the critical temperature T c . Second, 
this amounts to using the zero-tempcrature propagators. 
First results indeed indicate that transversal and longi- 
tudinal gluon and ghost propagators are little modified 



43l . |44j, |45[ for higher Matsubara frequencies 2itTn for 



n > 2,3. The biggest change appears in the gluon prop- 
agator longitudinal with respect to the heat bath that 
develops some enhancement compared to the transversal 
counterpart. The inclusion of the full temperature de- 
pendence is necessary for an accurate determination of, 
e.g., the critical exponents or the equation of state (see, 
e.g., [HI). This will be subject of a forthcoming paper. 



A SIMPLE ORDER-DISORDER CONFINEMENT 
CRITERION 

The preceding analysis gives rise to a simple confine- 
ment criterion which relates the IR behaviour of gluon 



and ghost 2-point functions to the deconfinement order 
parameter. Integrating the flow (|4|), we obtain 



T[A] = iTrlnr( 2 <°)[0,,4] +0{d t T%' Q) ) 



c.t. 



(14) 



where the counterterms (c.t.) denote the appropriate UV 
initial conditions of the flow, and the 0(8^^'°^) terms 
correspond to integrated RG improvement terms. The 
first term is explicitly regulator- independent, and so is 
the improvement term. This can be used to show within 
the specific choice © that the improvement term is sub- 
dominant for the following analytic argument, which is 
confirmed by the full numerical solution. 

The effective action in (|14p involves the Lapacian — D 2 
for vanishing field strength. In the constant Aq back- 
ground, we use the paramctrisation <?j4q = 27rT0 a , where 
4> a is a vector in the Cartan subalgcbra. The spectrum 
of the Laplacian then reads 

spec{-£> 2 L4 ]} = f + (27rT) 2 (n - ^|0|) 2 , (15) 

where the vg denote the N 2 — 1 eigenvalues of the her- 
mitian color matrix T a cj> a {T a ) bc = ~if abc being the 
generators of the adjoint representation. From Eq. (|15[) . 
it is clear that <f> is a compact variable. 

At high temperature, 2irT S> Aqcd, the effective po- 
tential is dominated by the perturbative regime, and the 
background-covariant inverse propagators of both gluons 
and ghosts are approximately given by their tree-level 
values r( 2 )- trco (-£> 2 ) = -D 2 . The perturbative limit of 
the effective potential V in d > 2 is given by the well- 
known Weiss potential |47[ , 



rVV fia 



1 1 

- + 2 



y ^Trln 



D 2 [A Q ]) (16) 



(d-2)r(d/2) 



T 



JV„ -1 
1=1 n=l 



E 



cos 2~Knvg\<p\ 



where the terms in curly brackets in the first row denote 
the contributions from transversal gluons, longitudinal 
gluons and ghosts, respectively. In the second row, we 
have dropped a T- and field- independent constant. The 
Weiss potential exhibits maxima at the center-symmetric 
points where L[(v4o)] = 0, implying that the perturbative 
ground state is not confining, (L) ^ 0. 

Now, wc perform the same analysis at low temperature 
ItxT <C Aqcd ■ The series in (|16|) converges rather rapidly 
due to the l/n d suppression of higher terms. Hence, 
the effective potential V{(p a ) is dominantly induced by 
fluctuations with momenta near the temperature scale 
p 2 ~ (2ttT) 2 . This does not change qualitatively in the 
presence of a non-trivial momentum dependence of the 
propagators. Wc conclude that only the first 10-20 Mat- 
subara frequencies play a role. Moreover, changing the 
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FIG. 2: Order-parameter potential for SU(2) (left panel) and SU(3) (right panel) for various temperatures. For SU(2) we show 
the potential for T = 260, 266, 270, 275, 285 MeV (from bottom to top). We find T c « 266 MeV for SU(2). In case of SU(3), 
the relevant minima occur in the Aq direction in the Cartan subalgebra. A slice of the potential in this direction is shown for 
T — 285, 289.5, 295, 300, 310 MeV (from bottom to top). A magnified view on the potential at the phase transition is shown 
in the inlay, revealing the lst-order nature of the phase transition with two equivalent minima at at T c m 289.5 MeV. 



propagator for the first two or three Matsubara frequen- 
cies, even though their weight is higher, only gives rise 
to minimal changes in the potential. This fully justifies 
the zero-temperature estimate on the propagators. 

With the parametrisation ©,([7]), the dressing func- 
tions Za(p 2 ), Zc(p 2 ) in the KOGZ scenario are charac- 
terised by the power-law behaviour (|SJ) in the deep IR, 
p 2 <C Aq CD . For low enough temperature, the spectral 
window — D 2 L4o] ~ (2ttT) 2 is in this asymptotic regime, 
and thus the effective potential arises dominantly from 
fluctuations in the deep IR, 

v m m = {^(1 + «a) + 5-(1 + «c)} 

xl T rln(-^ 2 [A ]) (17) 

= | 1+ (^-iW- 2Kc W uv (^). 

If the anomalous dimensions are such that the expression 
in curly brackets becomes negative, the effective poten- 
tial is reversed and the confining center-symmetric points 
become order-parameter minima. 

We conclude that the effective action (|T7)) predicts a 
center-symmetric quark-confining ground state if 

f{ KA , K c ;d) = d~2+(d~ 1)ka - 2k c < 0. (18) 

Provided that the 0(d t T k (2, 0)) terms in Eq. ([14]) remain 
subdominant, this equation provides a simple, necessary 
and sufficient criterion for quark confinement in Yang- 
Mills theory: if Eq. (fT8|) is satisfied the order parameter 
for quark confinement vanishes, (LL4o]) = 0. It is satis- 
fied for the whole one-parameter family of infrared solu- 
tions of Landau-gauge Yang-Mills theory. For the scaling 
solution with the sum rule ([9|), we are led to 

d — 3 ,„„. 

K=K C >^— ■ (19) 



which is satisfied for the numerical values for the scaling 
exponents k,j in rf = 2,3,4, see 0, HI). Specifically in 
d = 4, we have Eq. (fTU|) , and hence 

/(-2k c ,k c ;4) = -2.76.... (20) 

For the decoupling solution ifTTj) , we are led to 

/(-l,0;d)=-l. (21) 

Both values imply confinement, and hence the whole one 
parameter family of solutions is confining. Note that this 
is to be expected as corresponding propagators can be 
obtained within lattice simulations with different gauge 
fixings. 

The above confinement criterion has to be compared 
to the Kugo-Ojima criterion for color confinement k > 
and the Zwanziger horizon condition for the ghost k > 
and for the gluon k > 1/2 in d = 4. The Kugo-Ojima cri- 
terion and the Zwanziger horizon condition are necessary 
but not sufficient for confinement. Indeed for <k < 1/4 
in four dimensions, we observe that the Kugo-Ojima cri- 
terion is satisfied but does not lead to confinement ac- 
cording to the present confinement criterion (fT9| . We 
would also like to emphasise that, in effective theories for 
QCD, Eq. (fT8|) only serves as a necessary condition. It 
only restricts the propagators, and other Green functions 
in effective theories might violate related constraints. 

RESULTS FOR THE PHASE TRANSITION 

In contradistinction to the simple confinement crite- 
rion put forward above, the physics of the confinement- 
deconfinement phase transition, e.g., the transition tem- 
perature and the order of the phase transition, is deter- 
mined by the dynamics of the system and not by its IR 
asymptotics. Indeed, we find that fluctuations in the non- 
perturbative mid-momentum regime induce the center- 
symmetric minimum of the Aq potential long before the 
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propagators acquire their deep IR scaling form (jSJ). As 
only the deep infrared is sensitive to the infrared bound- 
ary condition the critical temperature is insensitive to 
this choice which is confirmed in the explicit computa- 
tion. 

The results presented below are achieved by numeri- 
cally integrating the flow equation (fj| in order to obtain 
the potential for an Aq background. The present trun- 
cation is optimised by using Landau-gauge propagators 
and RG improvement terms at zero temperature com- 
puted from the FRG for different infrared boundary con- 
ditions. It is also compared to results obtained by using 
fits to Landau-gauge propagators as measured by lat- 
tice gauge theory [7j and the RG improvement computed 
in [8j. For our numerical study of the order-parameter 
potential we have suitably amended the lattice propaga- 
tors by the perturbative behaviour in the UV and the 
corresponding power laws dHJ) in the IR. In Fig. Q] we 
show the gluon and ghost propagators as obtained from 
FRG computations || and lattice simulations ■ There 
is an impressive agreement of the results for the ghost 
and gluon propagators for momenta larger than about 
P ^ 700 MeV which holds for the whole one parameter 
family of solutions including the scaling one. The results 
for the ghost dressing from scaling solution of the FRG 
and lattice simulations start deviating for p < 700 MeV 
whereas the scaling solution for the gluon starts deviating 
for even lower momenta. Since the lowest non-vanishing 
Matsubara mode is associated with momenta at about 
\p\ ~ 2ttT c ~ 1700 MeV, the differences in the IR arc 
hardly probed in the present study of the deconfincment 
phase transition. This is confirmed by the explicit com- 
putation. In the vacuum limit, T — > 0, the picture arising 
from the preceding simple confinement criterion is con- 
firmed: a sufficient amount of gluon screening with or 
without an IR enhancement of the ghost creates a center- 
disordered ground state with quark confinement. 

The confinemcnt-dcconfincment transition is taking 
place in the mid-momentum regime that interpolates be- 
tween the perturbative regime and the IR asymptotics. 
The effective potentials for SU(2) and SU(3) for vari- 
ous temperature values near the phase transition are dis- 
played in Fig. [2l For SU(3) (right panel), the slice of the 
potential in A$ direction is depicted where the relevant 
minima for the phase transition occur. Reading off (Aq) 
from the minimum of the potential at a given tempera- 
ture, we can determine L[(A )] which is plotted in Fig. [3] 
For SU(2) (blue/dashed line), the phase transition is of 
second order. For SU(3) (black/solid line), we clearly ob- 
serve a first-order phase transition at a critical temper- 
ature of T c ~ 284 ± lOMeV with a lattice string tension 
^ = 440MeV, that is T c / ^ = 0.646 ± 0.023. The er- 
ror relates to the uncertainties of the fits for the lattice 
propagators which exceed the estimate on the system- 
atic error in the FRG computation. The result compares 
favourably both qualitatively and quantitatively with lat- 




FIG. 3: Polyakov loop for the Aq expectation value L[(3{Ao)] 
for SU(2) (blue/dashed line) and SU(3) (black/solid line). 
The phase transition is of second order for SU(2) and of first 
order for SU(3). 



tice simulations, see e.g. 0, EU. Also, our result for 
L[(Ao)] in the deconfincd phase is higher than the lattice 
measurement of the Polyakov-loop expectation value (L) 
in agreement with the Jensen inequality L[(^4q)] > (L). 
Note however that this statement has to be taken with 
care as the lattice result involves a non-trivial renormali- 
sation factor which is absent in the definition of L[(Aq)]. 
Indeed, L[(Aq)] < 1 whereas the rcnormaliscd Polyakov 
loop (I/) r en necessarily exceeds unity for some tempera- 
ture range as can be deduced from perturbation theory. 

As discussed above, corrections to our estimate arise 
from finite-T modifications of the propagators as well as 
from order-parameter fluctuations; the latter are more 
pronounced for SU(2) owing to the second-order nature 
of the transition. As expected, the critical temperature is 
not sensitive to the one-parameter family of solutions, it 
is only sensitive to the mid-momentum regime at about 
I GeV. Indeed, this also explains the fact that the gluon 
mass parameter is restricted from below: small gluon 
mass parameters would also trigger changes in the mid- 
momentum regime and almost certainly change physical 
quantities such as the critical temperature. 

In summary, we have established a simple confinement 
criterion that relates quark confinement to the infrared 
behaviour of ghost and gluon Green functions. This con- 
finement criterion is applicable in arbitrary gauges. Our 
full numerical analysis of the IR dynamics predicts a 
second-order phase transition for SU(2) and a first-order 
phase transition for SU(3), the critical temperature of 
which is in quantitative agreement with lattice results. 
The related Polyakov loop potential also plays an im- 
portant role for full QCD computations with dynamical 
quarks within functional methods, for first results on the 
QCD phase diagram see jiij . 
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